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ELF.CTROCHHMICAL MODELS FOR THE DISCHARGE CHARACTERISTICS 
OF THE NICKEL-CADMIUM CELL 

INTRODUCTION 


Basic electrochemical principles were applied to the sealed nickel-cadmium cell in an effort 
to describe its behavior and to predict its operating characteristics. An essentially thermodynamic 
approach was taken to arrive at several related but different equations describing its discharge. 

This investigation represents but one phase in an extensive effort to cliaructerize the 
behavior of the nickel-cadmium cell. The study presented here does not incorporate such opera- 
tional variablei as depth of discharge, long-tenn aging, n.*reatod cycling, preconditioning, rate of 
charge, etc. but merely considers the potent**?! Mine characteristics of a “preconditioned,” fully 
charged cel* dixharged at constant cunen' (. .«nsfant current operation was selected because it is 
simpler to describe and to facilitate compajison with the long-tenn experimental data which have 
been accumulated in that mode (1). 


MODEL I. SIMPLE THERMODYNAMIC MODEL 


The reaction taking place at the positive electrode of a nickel-cadmium e'en, i.e., the cathode 
during discharge, is usually taken to be (2.3): 


NiOOH + H,0 + e ^ NKOH)^ + OH~ 
Thc Nernst equation corresponding to this reaction is 
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where R is the gas law constant. F the faraday. n the number of electrons, and T the absolute 
temperature. Upon reananging terms. 
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At the negative electrode the reaction is (2) 


Cd(OH), + 2e CV + 2 OH" 


(3) 
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The corresponding Nemst equation is 
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or, upon rearranging 


E...„ = E” „ + — Cn 
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The Nemst equation for the net cell reaction would then be given by E^.gn = Ep^^ - or 

RT “NiOOH RT “Cd(OH)2 rj 


(4b) 
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Now let us make the following assumptions: 


“F (5) 

“Cd ^ 2 


(a) The Cd and the CdfOHlj at the negative electrode appear as separate solid phases 
so that a(^j = “ * 


(b) The NiOOH and Ni(OH)-, at the positive electrode are present in a single phase, i.e., 
a solid solution. 

Barnard, Randell and Tye (4,5) have presented evidence in support of this second hypothesis. Under 
these conditions the activities of NiOOH and NKOH)^ would not be unity but could be expressed 
as their respective mole fractions; 


“NiOOH “ ^NiOOH “Ni(OH)2 " ^Ni(OH)2 ' 
If NiOOH and Wi(OH)-> are the only species present in the solid solution, then 


^NiOOH ^Ni(OH)2 ~ ’ 

Representing the mole fraction of NKOH)t by X then the mole fraction of NiOOH would be I - X. 

The quantity X may also represent the state of discharge, e.g., X = 0 corresponds to the fully 
charged material and X = 1 to the fully discharged. In addition, if one combines the term involving 
activity of water and the standard potentials into a formal cell potential, E® , then the Nemst 
equation has the forni 


^cell 




(6a) 


Under conditions of cur.vnt How the cell resistance must be taken into account so that the 
discharge curve will have the fonn 


F 




(6b) 


where i is the total current, is the cell resistance and X is the fraction of discharge. 

The fonnal potential for the C'd(OM)-)/Cd system is reported (2) as - 0.869 volts vs. Hg/HgO. 
and that for the Ni(X)H/Ni(OII)% system as 0.365 volts vs. Hg/llgC) (2). Thus the formal potential 
for the net cell reaction is 0 365 + 0.864 = l .264 volts, and the discharge equation liecomes 


F = 1.264 + ^ Cn (7) 

A calculated curve based on equation 7 is shown in Figure I. An expenmental curve 
obtained during a capacity test in the (Vane series of nickel-cadmium battery tests (1) is shown for 
comparison. The cell resistance was calculated from the experimental potential at .-iOf discharge. 
Agreement with the expenmental data is not ver\- good. Further discussion of the errors will l*e 
pri'sented later. 


MOni L II. MOniFlFDTIlFRMOnYNAMICMOnFL 


Barnard, Randell and Tye (5) have descrilvd nickel hydroxide electrodes, based on evalua- 
*ion of the free energy of liiixlng two species in solid solution. The total free energy of the system 
is expressed as the sum of the free energy of reaction and the free energy ol mixing. In addition, 
a stvcalled “excess energy" tenn is included to account for interactions Ix'tween the species in the 
solid. M a result jf their treatment, they arrived at an expression for the potential of the nickel 
hydroxide system: 



A value of A/RT = 0 r\*presents "ideal” behavior and, in fact, corresponds to the simple themuv 
dynamic approach of Model I. N'alues of A RT = I or -I represent respectively positive and nega- 
tive deviation from ideality. The case where .A RT - 2 corresponds to the borderline kVtween a 
single-phase and a two-phase system. The reader is urged to consult the original paper lor details. 

When this modified thermodynamic approach is included, the equation for a NK'd dis- 
charge curve becomes 


i: = 


RT 

1.264 + — (n 
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Figure 1. Calculated discharge curve based on the simrl*.’ thermodynamic model, equation (7). 

The broken curve is experimental data. 

where K = A RT. The value of K was determined by a least s(iuares fit using four different dis- 
charge cimes obtained from the Crane data. The mean “best fit” value for K was 0.789. Thus, 
the discharge equation Ivecomes 

RT / 1 X\ RT 

I- = 1.294 + -^ Cn + 0.789 — f2X-l) - i Rq (9) 

.A calculated discharge curve based on equation 9 is shown in Figure 2. As before, an experi- 
mental curve IS included for comparison. It is immediately apparent that the calculated curve is in 
much closer agreement with the experimental curve than is that of Model I. The effect of diffeamt 
values ot current and temperature using equation 9 is shown in F'igua* .3. 
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l-igurc 2. Calculated discharge curve based on the modified thermodynamic model, equation (9). 

The broken curve is experimental data. 

MOni 1. III. MODII IH) THERMODYNAMIC MODEL WITH VARYING RESISTANCE 


In the previous model it was assumed that the cell resistance remains constant throughout 
the discharge. However, it has been shown that Ni(OH)'t in the pure state is a poorly conducting 

p-type semiconductor (6-9). On charging, n-type conductivity can be developed (8). It seems 
reasonable that the cell resistance would vary with the extent of discharge, reaching the lowest 
resistance when fully charged and the highest resistance when fully discharged. To this end, it 

was assumed a priori that the resistance takes the form R = Roe'^ * where X is the fraction of 
discharge. Thus R^, represents the cell resistance when fully discfiargcd. The disenarge equation 
then becomes 


CELL VOLTAGE CELL VOLTAGE 







0 10 :0 .<0 40 50 f>0 70 80 *^0 100 
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Figure 4. Calculated discharge curve based on Model 111. equation (10). 

The broken curve is experimental data. 


P = 1.294 + ^ Cn ^ ^ ^ (-X-I) - i Ry Kxp(-(l-X) 


As K'fore. K was treated as an adjustable parameter and evaluated by a least squares fit. The mean 
“best fit” K = The calculated discharge curve based on equation 10 i$ shown in Figure 4. 

As before, an experimental curve is included for comparison. 
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MODEL IV. CHRONOPOTENTIOMFTRIC MODEL 


“Chronopotentiometry" refers to an experimental technique in which electrolysis is carried 
out at controlled current and the potential it monitored as a function of time (10). The simples, 
and most common approach involves electrolysis at constant current. Under these conditions, the 
electrolysis reaction proceeds at a constant rate. If the reaction at one of the electrodes is the 
limiting factor, i.e. if one electrode is much smaller or if it has much less electroactive material 
available, then the shape of the potential-time curve will depend primarily on the reactions at that 
electrode. The potential assumes values characteristic of the redox couple and varies with time as 
the oxidized/reduced concentration ratio changes at the electrode surface. Eventually, after the 
concentration of reactant drops to zero at the electrode surface its flux is insufficient to consume 
all of the electrons which cross the electrode-solution interface. In the case of a reduction, the 
potential will then rapidly shift toward more negative values until a new second reduction process 
can start. The time required for this potential transition to occur is called r, the transition time. 

A typical chronopotentiogram is shown in Figure 5. 



I'lKurc 5 A typicul chronopolentioprani. file quarter-wave potential. 
l>/4 . correspv>nils to tlic potential .it a limc » * I '4 r. 


K 



If ebctron transfer is rapid, the Nemst equation applies. It can be shown (II) that the 
potential*' ime curve has the form 


E - E^/4 



Rearranging the terms within the logarithm one obtaiiu 


E “ Ef/4 


RT 
+ — 
nF 


fin 



(lie) 


(llh) 


Since electrolysis is being carried out at constant current the amount of charge pas^d at nny time, 
t, after the start of electrolysis is given by the current time product, q « it. In addition, th*^ chaige 
corresponding to the total capacity of the electrode is given by the current and the transition time, 
Q * it. Thus, the fraction discharged X ■ q/Q * t/r. In addition, % E°' (II). 1 .icn the 
lischarge equation is 


E 



( 12 ) 


A calculated discharge curve based on equation 1 2 is shown in Figure 6. An experimental curve is 
included for compar son. 



Figure 6. Calculated discharge curve based on Model IV, equation (12). 
The broken curve is experimental data. 
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MODEL V. MODIFIED CHRONOPOTENTIOMETRIC MODEL 


In the normal chronopotentiometry experiment the current is held constant thus forcing 
reaction to occur at a fixed rate. The presence of a finite capacitance, which might arise from the 
double layer at the electrode^Iution interface or within the deposit of active material, introduces 
a non-faradaic charging current, i^, proportional to dE/dt. Thus, only a portion, if, of the total 

current goes to the faradaic reaction 


i = if 




(13) 


Even though i is constant, and if vary with time. This effect leads lo a distortion of the potential- 
time curve, de Vries (12) has shown that the charging contribution is represented by the dimen- 
sionless parameter 

K - Cd 

nFC(,(»[)T)^ 

where Cj is the capacitance, C* the bulk concentration and D the diffusion coefficient. Also 



n F A d'^ 


* 



0 


which is known as the Sand equation (13). Rearranging equation 14a one obtains 


nFD'4w»''r* = — 
^ nF 



(14b) 


(14c) 


and from equation 14b 



n F C, 
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Then, 


- 'f _ 1^ 
A nF 


(14d) 
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Solving for 1 It, 


(140 


1 

T 



‘f 


and 


where X = ift is the fraction of discharge. Substituting equation 1 S into equation 1 2 and letting 
(K/CjjA) = K one obtains 



(16) 


The value of K was evaluated by a least squares fit of the experimental discharge curves. The mean 
best fit value for K was 0.01 208. The equation for the discharge curve then becomes 


RT 

E = 1.294 + — in 
F 


1 - (0.02416— X) 
\ RT / 


>/2 


- i Rq 


The calculated discharge curve based on equation 1 7 is shown in Figure 7. 



*1 mtrHA»( r 

Figure 7. Calculated discharge curve based on Model V, equation (17). 
The broken curve is experimental data. 


(17) 


COMPARISON OF THE MODELS 


The equations corresponding to the models discussed in this report are summarized in 
Table 1. 

For each model, an error distribution was calculated with respect to each of the four experi- 
mental discharge curves. A typical set of error curves is shown in Figure 8. A summary of the 
respective standard deviations is shown in Table 2. As might be expected, Model I is in poor agree- 
ment with the experimental data (see Figure 1, Figure 8 and Table 2). None of the calculated 
discharge curves is in complete agreement with experiment for the full discharge. For Models II, 

III IV and V, agreement seems quite good for 50-100% discharge (see Figures 2, 4-8 and Table 2). 
Standard deviations calculated for the full curve, 0-50% discharge and 50-100% discharge are 
included in Table 2. In each case, however, agreement is poorest in the initial portion of the dis- 
charge curve. 


Table 1. Summary of Equations for the Five Models 


Model 

1 Simple Thennodynamic 

Equation^ 

RT /1-X\ 

E=,.294 .-«„ (— )-iRo 

11 Modified Thermodynamic 

RT 

E = 1.294 + Cn 1 

nF 

^ l-x\ RT 

. 0.789 -(2X-1) 

- i Rq 

ill Modified Thermodynamic 
with Varying Resistance 

RT 

E = 1.294 + — in 
nF 

( 1-X\ RT 

- iRfle-d-X) 

IV Chronopotentiometric 

RT i 

E = 1.294 + Cn| 

nF ' 


V Modified 

Chronopotentiometric 

RT 

E = 1.294 + -— Cn 
nF 

1 -(.01208 |txf 

- i Rq 

( 2nF 

(.01208 X ) 

V RT ' _ 


is the fraction of discharge, Rq is the cell resistance 


1 
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Figure 8. Typical error curves for the five models. 

The curves shown are all with respect to the same experimental discharge. 
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Table 2. Comparison of Errors for the Models 


Standard Deviation^, millivolts 


Model® 

0-100% 

0-50% 

50-100% 

1 

14.4 

13.9 

14.8 

11 

8.1 

9.5 

5.7 

III 

9.7 

11.1 

7.6 

IV 

13.2 

14.0 

12.0 

V 

10.8 

14.0 

6.4 


^Correspond to the equations listed in Table 1. 

^Each standard deviation is the mean for the four experimental discharge curves 


None of the models accurately represent the initial portion of experimental discharge 
curves. The disagreement may be the result of an artifact in the experimental procedure used to 
collect the discharge data or perhaps the models neglect some fundamental process which occurs in 
tlie early stages of discharge. Tliis problem certainly warrants further investigation. 

Models II and 111 appear to agree more closely with experiment than do the others. Best 
agreement in all regions of the discharge curve is found with Model II. 
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